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1. Introduction. 

The problem of estimating energy exchanges between fast and slow degrees of freedom is a very relevant 
one in many domains of physics, and provides a main motivation for studying adiabatic invariants. Very 
elementary models come from gas dynamics, where the slow degrees of freedom correspond to the motion 
of the molecules center of mass, and the interaction with the fast internal degrees of freedom is provided 
by molecular collisions. In the simplest case one has just one vibrational and one translational degree of 
freedom; the Hamiltonian is: 

H{x,y,p,,Py) = ^ + Voix) + ^{^+^iu^y^) + Viix,y) , (1.1) 

and represents the collision of a diatomic molecule with a fixed wall, or with a point mass, in one dimension. 
Here x and are the canonical coordinates of the center of mass, while y and py are the coordinates of the 
internal degree of freedom; m is the total mass, /j, the reduced mass and w/27r is the vibration frequency. 
Both Vo and Vi are supposed to be smooth (in fact analytic) and to vanish for \x\ co. One can take 
Vi(x,0) = 0; at fixed total energy one has then \y\ ~ and thus Vi is small for large u), (generically) 

Vi ~ L0-\ 

Models like (1.1) have been introduced by Jeans [J1,J2] to investigate the times of relaxation to equilibrium 
in polyatomic gases, and to understand why in ordinary conditions, and in apparent confiict with the 
equipartition principle, the vibrational degrees of freedom of diatomic gases do not take part to the energy 
exchanges. In fact it is not well known that Jeans conjectured (after Boltzmann [Bol,Bo2]) that the apparent 
"freezing" of fast vibrations, in this and other problems, is simply due to the fact that the energy exchanges 
with fast degrees of freedom are too small to produce any appreciable effect in a reasonable time scale (see 
[Bel,Be2,BGGl] for comments on this point). On the basis of heuristic arguments, Jeans proposed for the 
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energy exchanges an exponential law of the form 

~ 5e-™ , (1.2) 

where £ and r are natural units of energy and time (of the order of the energy per degree of freedom and 
of the collision time, respectively). As far as we know, this is the earliest proposal of an exponential law in 
connection with a problem of adiabatic invariance. Later Landau and Teller [LT] reconsidered the question 
in connection with the problem of sound dispersion, and worked out an elementary approximation scheme, 
actually not so different from Jeans' original one, which, once applied to models of molecular dynamics like the 
above (1.1), see [R], immediately leads to exponential laws like (1.2). We shall refer to this approach as to the 
Jeans-Landau- Teller (JLT) approximation. Besides the original papers [J1,J2,LT], the JLT approximation 
is described in [R,BCS,C]. 

To recall the essential result, following [ECS], we use for the fast degree of freedom the action-angle 
coordinates I, Lp, and consider a general Hamiltonian like: 

p 

Hix,p,if,I) = h{I) + — + Vo{x,r,e) + eVi{x,I,if;e), p = Px ■ (1-3) 

For a vibrational degree of freedom, which is the case in which we are mainly interested, one has h{T) = loI 
and (generically) e ~ bj~^ . The coupling term Vi can be conveniently taken to have zero average in ip, as 
the average could be added to Vq. 

A possible formulation of the JLT approximation is as follows: 

2 

i) Ignore Vi, and solve the equations of motion for the reduced Hamiltonian Ho{x,p; I) = ^ -|- Vo{x,I), 
where 7 is a parameter; let x(t), p(t) be the solution for given initial data (a;°,p°) at t = 0. 

ii) Use this solution, together with the free motion I(t) = 1° , (f(t) = ip° -\- uit, where ui = ■|j(7°), in the 
Hamilton equation for 1:1= —e^^(x(t), I{t), (pit); e), integrate over time and compute Ai? = h{I{-\-(X))) — 
M/(-oo)). 

As a result (see [BCS] for details), one gets AE" in the form of a Fourier series in Lp°: 

AE = e''"^" (1.4) 

with separate exponential laws for the different Fourier components: namely^) 

1^,1 = 5,e-l''l™ fori^^O, |^o| = ^oe-2™ (1.5) 

(proportionality to a power of w is possibly hidden in the constants the crucial constant r is given by 
the width of the analyticity strip of the Fourier coefficients of Vi(x(t), I(t), cp; e), thought of as functions of 
complex time t. In typical cases, and in particular for the potentials used in [LT,R], r is the width of the 
analyticity strip of the approximate solution x(t), p(t), and is easily computed. The precise identification 
of r, as well as the result that the average energy exchange Eg decays much faster than the zero-average 
part (the latter practically coincides with the first Fourier term), are the essential achievements of the JLT 
approximation. 

As shown in [BCS], the JLT approximation, although apparently crude, works indeed excellently, and in 
particular the separate exponential laws for the different Fourier components, as well as the identification 
of r, are fully confirmed by accurate numerical computations (the agreement between the theoretical value 
of r and the numerically computed quantity is better than one percent, for relative energy exchanges in 
the rather wide range 10"'* ^ 10"^*^). The reason of such a good agreement is not so clear; in particular 
(see [CBG] and [BCS] for comments) one does not understand why the analyticity strip r of the zero-order 
solution x(t), p(t) plays such a relevant role. One should note that, apart from exceptional cases, r is not 
an approximation to the analyticity strip of the true solution: the latter is totally different from r, and even 
vanishes at large uj [BCS]. 

The purpose of this paper is to understand a simplified version of the problem (1.3) from a mathematically 
rigorous point of view. We shall restrict ourselves to the case of a harmonic oscillator, namely h(I) = ujI 

Properly speaking, in the exactly isochronous case h{I) = a;/, at first order in e the average Eq vanishes: contributions 
proportional to e~^'^^ are however expected at higher orders. 
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and to /-independent Vi, Vo, and work perturbatively with e as the perturbative parameter; only at the end 
we set s ~ u)~^. Although our approach is in principle more general, we explicitly consider here the case 
Vo{x) proportional to e~^/' where / is a length scale; moreover, we shall assume that Vi(x,Lp,e) factorizes 
as g(ip)Vi(x), with Vi(x) decaying sufficiently fast at cxd; we shall also occasionally refer, for definiteness, to 
Vi proportional to e~^/'. More general choices of Vo and Vi could be allowed at the only expense of a more 
involved notation: some comments on such generalizations will be made later, but we prefer to sacrifice the 
generality of the model to the exposition of a general technique in a rather special class of problems. We 
have not attempted, however, to replace the coupling strength e in (1.3) with a power of ^/I and to take 
Vo, Vi analytic in ^/I (rather than I): the latter problem may require new ideas. 

Let us stress here that, although we treat explicitly the case of a single fast variable, our approach can be 
generalized to the case of several fast variables; for a preliminary result in this direction see Section 4.4. 

Concerning the technique, we shall proceed as in [CG,G3,G1,CF1], and expand solutions in series of e, by 
using a diagrammatic technique to represent AE and to look for relevant compensations. The first order 
of this expansion gives the JLT approximation, and due to the decay property of Vo and Vi at infinity, the 
perturbative series converge. In particular we obtain the Fourier coefficients Ei, of the energy exchange in 
the form of a series in s: E^, = X]ft>i ^^^^h,!/, and (as is very crucial) all terms are bounded by the same 
exponential factor: precisely we get 

Eh,, S-'^'+^e-^"^'^^^^-^^ for 7^0, Eh,o 6-'"'+^ e-^'^^^^-^\ (1.6) 

with some given q > depending on the model, while (5 > is a small free parameter. The leading terms 
Ei,±i are written in integral form. Denoting by Vi_y the Fourier components of Vi it is 

/oo 
Vi,±i{x{t),T{t))e^'^'dt (1.7) 
- OO 

In some interesting cases such terms can be explicitly evaluated: for example, for the very special choice (see 

Vo(a;) ~ e"^/' , Vi(a;, ^, 7) ~ e"^/' cos ^ , />0 (1.8) 

and in this case it is Ei,±i ~ u)^e~'^'^ . For such a model one finds q = A and r = 'iT/(2^/rf), rj" denoting the 
asymptotic value of the translational energy at time — cxd. 

Inequalities (1.6) clearly solve Jeans' problem of giving an exponential upper bound on AE. Indeed, taking 
for example s = I/toj and making the choice b = e'^ with Q < ^ < ,hy summing the series one gets that 
the first term dominates (for the values of Lp for which it is not zero) for small e and 

\AE\ < const. , = r[l - (rc.)-^] , (1.9) 

(TUl ) " 

SO the equilibrium times are necessarily long (see [S] for preliminary results and comments on this point). 
Unfortunately, from the above analysis one cannot conclude that the first order in e, namely the JLT 
approximation, dominates the expansion. Hence we cannot deduce from (1.6) a two-sided (or "exact") 
estimate for AE: in fact the correction I — 6 at the exponent implies that the other (higher order) terms 
could exactly cancel the first order contribution, leading to a much smaller AE; or they could even dominate 
leading to a AE much larger than the one of the JLT approximation (but of the order of 0(e~'^^ '••^~^^), 
of course). A two-sided estimate can be obtained if one assumes £ ^ in this case one can then take 

6 ~ and the correction 1 — (5 is then rigorously negligible. 

Remark: In special cases, for instance with tidal dipolar forces in (1.1) (overall potential Vo(x) + Vi(x, y) = 
V{x + ij/) + V{x — \y), with some V), one has Vi{x, y) = 0{y''^), and thus s ~ but this is not enough 
as one should have Vi = 0{y^'^) in order to be able to justify a choice of e of order 0(w~*) to be able find 
upper and lower bounds on AE with the above estimates (1.6). The only case in which (1.6) suffices for 
getting sharp upper and lower bounds seems to be the case in which the amplitude \fTJu) of the oscillations 
is initially sufficiently small (one could see that y remains then consistently small at any time). In such a 
situation it is reasonable to assume e to be small compared to and deduce from (1.6) a lower estimate, 
too. This can be interpreted as indicating that an initially unexcited molecule may reach an energy of the 
order within a computable, exponentially large, number of collisions, {i.e. an exponentially large time). 

A proof of this fact is, however, not a well posed question in the frame of the above models because it requires 
studying models in which several collisions can be observed (while model (1.3) only considers one collision). 
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Our perturbative approach to the problem is quite different from the Nekhoroshev-like approaches to 
similar problems, see [N1,N2,BGG2,BF], which also lead to exponential estimates for the energy exchanges. 
Here we do not introduce chains of canonical transformations, rather (as in Lindstedt method) we expand 
solutions. The essential quantity, namely the exponential e~'^^ , appears already at first order in e, and not 
after optimizing a divergent series truncation. Moreover, the crucial constant r is precisely identified, and 
associated with an easily computable parameter of the unperturbed Hamiltonian; the corresponding quantity, 
in the Nekhoroshev-like approach, appears instead indirectly, through the analysis of an apparently divergent 
rate of perturbative series, and a sharp determination of it seems practically impossible. We also remark 
that here the results are more detailed (one explains in particular the different exponentials in the bounds 
for the average and for the different Fourier components). 

Our approach is instead rather close to the modern methods of investigating the problems of homoclinic 
splitting in forced pendula, and similar problems which are the basic steps to understand the Arnol'd 
diffusion, see [ACKR] and the review paper [Gl]; some comments on this point are deferred to Section 
4.3. A connection also exists with other versions of perturbation theory "beyond all orders", like the field 
theoretic approaches, see [PV,G2], and the modern KAM theory as started by Eliasson, see [E,G3,CF1]. 
Finally, our paper is quite close to [C], where one considers the JLT approximation scheme, and deduces an 
upper bound by a different technique, which essentially consists in a rigorous version of Jeans' original ideas. 

For a study of a quantum problem related to the problem discussed here, see [JKP], [JPl], [JP2], [MN]. 

This paper is organized as follows: in Section 2 we introduce precise notations, and state formally our 
results, limiting ourselves to the restricted case described above. Section 3 is devoted to the proof, which 
essentially follows [G1,G3]; the only new ideas are the use of the energy-time coordinates, in place of x,p, 
for the slow degree of freedom, and the particular "cancellation mechanism" which at any order h leads 
to the exponential bound (1.6) for Eo h- The latter is the only pomt where the Hamiltonian character of 
the problem is used in an essential way. Section 4 shortly accounts for some possible generalizations of 
our approach: precisely, in Section 4.1 we extend the results to a wider class of perturbations; in Section 
4.2 we show that, to get cancellations, the symplectic symmetry of the equations can be replaced by other 
symmetries, like the time reversal invariance of the equations; in Section 4.3 we establish the connection 
with the problem of separatrices splitting; in Section 4.4 we briefiy treat the case of several fast variables; 
and in Section 4.5 we shortly indicate a possible way to improve our estimates in some special cases. Finally, 
in the Appendix we illustrate how one could proceed without the use of the energy-time variables. 



2. Statement of results. 

To stress the mathematical character of our analysis we assume that all physical quantities are expressed 
in natural dimensional units; in other words, from now on we deal with dimensionless quantities. 
We consider the Hamiltonian: 

2 

H(x,p,^,I) = L,I+^ + Vo(x) + eg(^)Vi(x) ; (2.1) 

although we have in mind the general problem of scattering by a smooth potential barrier, we explicitly 
treat the case of Vq as in (1.8): Vo(x) = . Concerning Vi, we only need it be real analytic and to decay 
to zero for x +oo; for definiteness we shall occasionally refer to the simple choice 

Vi{x) = e-'' (2.2) 

(such exponential potentials are taken from [LT,R]). 
The one degree of freedom Hamiltonian 

2 

Ho(x,p)=^ + Vo(x) , (2.3) 

giving the motion of the center of mass for £ = 0, only admits unbounded (i.e. scattering) solutions, namely 
solutions with x(t) co and p(t) ipoo for t ±cxd for all energy values. We therefore introduce the 
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energy-time variables of the Hamiltonian (2.3), which are the analog of the common action-angle coordinates 
for unbounded motions. To this end, let us fix arbitrarily an origin x° for x, and consider motions with 
initial data x° and p° = \/'i{r\ — Vo(a;°)), namely motions with total energy r\. The corresponding solutions 
can be expressed as functions of r\ and of time: which we shall denote by ^. We write such motions as: 

X = X{^,T]) , p=p{^,T]) . (2.4) 

For the above choice of Vq one easily finds explicit expressions in terms of elementary functions, namely 

x(i,rj) = log[^{cosh^2(i-C)f] , p(i,rj) = ^ts^nh^2(i-C) (2.5) 

(with ^° such that x(0, rj) = x°). In general, in any scattering problem, the solution (2.4) is implicitly defined 
just by a quadrature and an inversion. But the analysis of the analyticity properties may be involved and 
difficult: and this is the main reason for restricting our analysis to Vq = e~^. 

Expressions (2.4) can be regarded as a change of variables from (x,p) to (^, rj), which is, of course, a canonical 
change of coordinates. After the change of coordinates (2.4), the Hamiltonian Hq becomes Ko(C,f]) = f], 
while the Hamiltonian (2.1) is changed into: 

K(C, r,,^,l) = ujl + r, + efii, r,)g{^) , (2.6) 

with 

m^) = Vi{x{^,^)) . (2.7) 

^From such Hamiltonian one gets Lp = to, and thus Lp{t) = Lp° -\- tot; for ^, rj one then finds the pair of 
time-dependent equations 

l + £/e(€,»?)ff(^''+^0 > r? = £/,(e,r?)ff(^''+c.t) , (2.8) 

with 

h = %^ h = -%- (2.9) 

The form (2.9) of expresses the Hamiltonian character of the problem. This however plays no role 

in the perturbative construction, and it is useful only at the very end, to see that many individually large 
additive contributions to the quantities that we want to estimate {i.e. energy exchanges between internal 
and external degrees of freedom) exactly compensate. We shall proceed with generic and , and recall 
(2.9) only when necessary. 

We now explain the assumptions we make on g, and . Concerning g, we assume it to be analytic in a 
strip jlmi^l < p, with some p > 0, and bounded by 

\g{^)\<l, for |Im^|<p (2.10) 

(taking one on the r.h.s. of the inequality is not restrictive). Properly speaking, in the special problem 
in (2.1) with perturbation independent of I, analyticity is not really necessary, and leads only to a minor 
improvement of the estimates. Analyticity is however necessary in the general case (see [N2] for a similar 
situation). 

Concerning and let us preliminarily consider the Hamiltonian case, with Vi given by (2.2). In such 
a case the function / defined by (2.7) is given by 

fit V) = ^ , (2.11) 

[coshy^(e -e)]2 

and for any real rj, f is analytic, as function of ^, in the complex strip |Im^| < r(ry) = By moving rj in 

a complex ball around a real value rj" , the singularities of / move, but / remains analytic in domains like 

Vsi^n = {(€, r?) e (T : |ry - ry"! < ; jlm^j < (1 - 6)t} (2.12) 

with r = T(ri°), 6 < I and a suitably chosen c. Simple considerations lead to estimates on /„, a = ^, rj: 

1 



j\(m - jy. 



^^lh_{s + ia,t)\<C^6-'^-'^"w{s) for \a\<{l-6)T{t) (2-13) 
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where mo = 3, C is a suitable constant, and w is a suitable integrable function: 

w{s)ds = W < oo (2.14) 



Under general assumptions on namely analyticity for ry in a small ball around a value rj" and ^ in a 

complex strip, and sufficiently fast decay (faster than l/|Re^|) for |Re^| ±cxd, one gets estimates of the 
form (2.13), (2.14), with suitable mo, C, W , and r = T(ryo)- 
We shall use mo, C , W and r as the quantities characterizing our problem, and express our results in terms 
of them. Let us note that, since ^ is the time of the unperturbed motion, the quantity r appearing in (2.13) 
is precisely the quantity appearing at exponent in the JLT approximation. 

For £ = the equations of motion (2.8) trivially admit the free solution S,o{t) = t, rio(t) = rj" . For e ^ we 
are interested in real solutions such that, asymptotically for t -co, 

C(t)-t^O, v(i)^v" (2.15) 
Therefore we introduce the Taylor expansions (which will be proven to be convergent) 

^(t)=t + J2^''Ui) , ry(t) = ry<' + ^£S,(t) , (2.16) 

h>l h>l 

and rewrite this, in compact notation, as 

a{t) = aoit) + ^e''ahit) , a = ^ or ry . (2.17) 

h>l 

The overall variation of a = ^ or ry (the asymptotic delay and respectively the energy exchange) is then 
= J2h>i £''aft(+oo). 

Denote by , i> ^ !E, the Fourier components of Uhi+co), a = ^ or ry, namely Uhi+co) = X^i^g^ s'"^ 
then the Fourier expansion of Aa is 

Aa = ^ A« e'''^° , with A« = ^£''aft,, . (2.18) 

v^'Sj h>l 

Let us remark that the asymptotic energy exchange AE" we discussed about in the Introduction, coincides 
with — Ary, and thus 

Eh,. = -ri^,., E, = -Al = -J2£''Vh,.- (2.19) 

h>l 

We shall prove the following 

Proposition 1. Consider the equations of motion (2.8): 1st p, mo, C, W and r be such that inequalities 
(2.10), (2.13) and (2.14) hold, for < 6 < I, and consider motions satisfying the asymptotic conditions (2.15). 
Then: 

i) For small e the Taylor series (2.17) for a(t) converge, precisely for any real t one has 

\an{t)\ < AB''-^ (2.20) 
with suitable constants A, B depending on C and W . 

ii) The Fourier components a^^v of ah{-\-oo) satisfy the exponential estimates 

\ah,„\<AB^-^8-^'i+^e-\''\'^''^^-'''^-p\''\, a = £.,r] (2.21) 

with q = mo + 1. 

Hi) In the Hamiltonian case (2.9), the average rjhfi of the energy exchange satisfies the exponential estimate 

\T]h,o\ < AB''-H-'"'+^e-^'^^^^-^^ (2.22) 
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iv) At first order in e one has the explicit expression 

/CO 
fa{t,V°)e'''^'dt , v^m. (2.23) 
- oo 

Possible values of the constants A and B are A = W , B = 8CW. 

Point iv) of the proposition reports the expression of the Fourier coefficients of A^, Arj at first order in e; 
as will be clear from the proof, one could produce explicit expressions for higher orders, too. For the choice 
(2.2) of Vi, one has in particular 

2 2 

nii=47riai , A£' = 87r sin^° + ... (2.24) 

with AE" determined up to corrections of order 0((£w*)e~^'^'^). 



3. Proof of Proposition 1. 
3.1. The recurrent scheme, and proof of item i). 

By introducing in the equations of motion (2.8) the Taylor expansions (2.16), collecting then terms of the 
same order h in e, one gets 

ah(t) = Fa^h(t) , a = (,r], h>l, (3.1) 

with 

F„,i(t) = g(^" +iot)U(t,rj") , (3.2) 

and for h > I 

l<m<h-l 0<j<m 

J2 ^kAt)--<kAi)m,+At)---mJt) 



\k\=h-l 



(3.3) 



where use has been made of Co(i) = i, and 

i!(m-i)! 5ei5r?™-J ' ' ' ^ ' ' 

Equation (3.1) is trivially solved by 

Mi) = aft(0)+ / Fa^h{t')dt' . 
■Jo 

The value of aft(O) is determined by imposing the asymptotic condition (2.15); one gets 

ah(t) = I F„^h(t')dt' . (3.4) 



Equations (3.2), (3.3) and (3.4) provide a recurrent scheme to compute C(t), ri(t) to all orders in e. 

We can now deduce the estimates (2.20), which in turn imply the convergence of the series (2.16). Using 
the assumption (2.13) for 6 = I, as well as |(j'(',c)| < 1, and denoting 

Ch = sup max {\ih{t)\, \Vh{t)\} , 
ten 
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we immediately get the recursive estimates 

Ci<W , Ch<W J2 (^ + 1)^^™ E Ck,--<k^; (3.5) 

\k\=h-l 

the estimates (2.20) are then an immediate consequence of the following lemma: 
Lemma 1. The numbers Zh, I < h < co, defined recursively by 

Zh = a E 6™ E Zfcj---Zfc^, zi = c (3.6) 

satisfy the inequality 



|fe|=^-i 



Zh < . (3.7) 

To prove (2.23) one applies this lemma with c = a = W and h = 2C; m + 1 in (3.5) is overestimated by 2™. 
One immediately gets (h < Zh < W(SCW)''-'^ /h, even stronger than (2.20). 

Lemma 1 is in turn an easy consequence of the following combinatorial Lemma 2, which also will be useful 
in the following. 

Lemma 2 The numbers Mh, I < h < co, defined recursively by 

Mh= J2 J2 Mk,---Mk^ , Ml = 1 , (3.8) 

l<m<ft-l fti,...,ftm>i 

\k\=h-l 

are given by 

2''-i(2/i - 3)!! 4''-! , , 

Mh = \, < ^ (3.9) 

Lemmas 1 and 2 are left to the reader with the hint that (3.8) implies for the generating function u(x) = 
X]ft>i MhX^ that it has to verify the relation u{xY = u{x) — x, m(0) = 0, which is solved by u(x) = 

^{^~-^/^-'ix). 

3.2. The tree expansion, and proof of points ii) and iv). 

Equation (3.4), with F^^ given by (3.3), gives aft(t) as a sum of integrals, one for each choice of m, j and 
ki, . . . , km, within the range specified in the corresponding sums. To each term one naturally associates an 
elementary tree (see figure 1), namely a root which ends in a vertex v and there divides into m branches; the 
root has a label t ^ IR and a label a = ^ or ry, while the branches carry labels a^^, . . . , a^^ (i-e., either ^k, or 
rjkj; the rule is that "^-type" branches must precede (top to bottom, say) "ry-type" branches. The labels a 
and t attached to the root, together with h = I + ki, indicate that the graph is a contribution to ah(t); 
the vertex represents an integration over a dummy variable ty, in the range — cxd < ty <t; the function to be 
integrated is the product of several factors, namely one factor Ck,(^v) or rj^^ity) for each outcoming branch, 
and an additional factor g{ip° + uity)f^'^ (ty , rj") for the vertex, j being the number of ^-type branches. For 
example, the elementary tree drawn in figure 1 represents a contribution to ^sii), namely 

g{^° +ioty)f^'^(ty)C2(tv)Ci(tv)v4(iv)dty . (3.10) 




Figure 1. An elementary tree contributing to 



29agosf 0I995 



8 



^From elementary trees one then constructs trees, just by recursively expanding all branches in elementary 
trees; the procedure stops when the end branches represent either or rji. The latter functions, in turn, 
are expressed via (3.4) and (3.2) as integrals of g and or Each tree clearly represents one term of the 
sum obtained by recursively expanding each of the ^k, or rjk^, appearing in (3.3), via (3.4) and (3.3) or (3.2), 
till everything is explicitly expressed in terms of g, and of the Taylor coefficients of . See figure 

2 for an example, actually a possible expansion of the elementary tree drawn in figure 1. Each tree 6 with h 
vertices (including the end vertices) has labels t G iR and a = ^ or ry on the root, indicating that it represents 
a contribution to ah(t); each branch in the tree has in turn a label a = ^ or ry, and the rule is that ^-type 
branches must precede ry-type branches issuing from the same vertex. The number of branches issuing from 
V is denoted niy, while the number of ^-type branches is denoted jy (jy = niy = for end vertices). 




Figure 2. A possible expansion of the elementary tree drawn in Figure 1. Labels of branches 
and numbering of vertices, corresponding to the multiple integral below, are reported in 
the drawing. 

The contribution, or value, of a tree is obtained as follows. Any vertex v ^ 6 represents an integration 
over a dummy variable ty, so for each tree with h vertices one has a multiple integral of multiplicity h. The 
quantity to be integrated is the product of one factor for each vertex v ^ 6, precisely 

gi^''+i,ty)fZ^-^^ity,7j'') (3.11) 

(g(ip° + i^ty)fay(iv , for end vertices), Uy being the label of the branch ending in the vertex v. Finally, 
the domain of integration respects the partial ordering^) of vertices in the tree: if t = (to, . . . ,th-i) G is 
the set of integration variables, then the domain is 

Ti(6) = {t = (to, ■ ■ ■ , th-i) G IR'^ '■ ty < tyi if v' precedes v ; to <t^ . 

For example, the contribution of the tree drawn in figure 2 is 

ft fto rti fto rto rti fte 

dto dti dt2 dta dt4 dt5 dte dtr g(^'' + cuto) ■ ■ ■ 

■CO J — oo J — oo J — oo J — oo J — oo J — oo J — oo 

■ ■■g(^" +iotr)f^''(to)fl°(ti)f,(t2)Mts)f^'HU)Mh)f/(t6)m7) 

(the way the integrals are nested refiects the tree structure). In fact, our aim is to compute ah(+(X)), and to 
this purpose we must extend the outermost integration to +cxd; this means that the domain of integration 
isToo(^?). 

We shall denote by 0^ the set of topologically distinguishable trees with h vertices (the order in which 
the branches appear from top to bottom is relevant, so the trees of figure 3 are different); the set of labels, 
namely a = {uy , v ^ 6} (the first component is the label a of the root) constitute the decoration of the 
tree. To prevent confusion, let us stress that a refers to the whole set of labels on the branches, including 
the root, while a refers to the root, and indicates the quantity one wants to compute. For each tree, the 
decoration a determines in particular j = {jy , v ^ 6}. Having introduced this formalism, we can write 



(3.12) 

V(e,a) = [ ■ ■ ■ [ T\g(^"+cuty)fZ^'^^(ty,^")dt; 



Any tree— like graph clearly induces a partial ordering in the set of vertices (from left to right in our figures). 
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the primed sum over a is intended to run over all allowed choices of the branch decoration, of course not 
on the label a of the root. The quantity V(6,a) represents the value of the decorated tree (0,a). Let us 
remark that the set of trees for and for rjh differ only by the label a = ^ or ry on the root. 




To separate the different Fourier components of Uhi+co), we must expand each function g appearing 
in (3.12) in Fourier series, and conveniently collect terms. To this end, we add further decorations to the 
tree, namely we attach to each vertex v ^ 9 a label ny G !E, corresponding to the chosen Fourier component 
of g, and collect decorated trees such that ny = i> (trees with momentum m, in the terminology of [Gl], 
[G2]). So we write 

with 

V{e,a,n)=(j{g^^) j ■■■ j [ H /™:''"(^- di (3.13) 
(the dot at the exponent denotes the usual Euclidean scalar product). Correspondingly one has 

"ft,- = E E' E • (3-14) 

We can now proceed to the estimate of a^^^, in order to prove point ii) of the proposition. To this 
end we consider the multiple integration on t in (3.13), and simultaneously raise all integration paths to 
Imty = ±r(l — 6), with sign independent of v and coinciding with the sign of wi^. The only "vertical" paths 
one introduces in this way are at infinity, so they give no contribution to the integral; the analyticity of all 
functions in the strip |Imt| < r, implies that the value of V(6,a,n) does not change. Having raised the 
integration paths, one immediately gets the exponential factors one is looking for, since one has 



e 



-\^\u,t(1-S) ^iwn^^ s = Ret 



The estimate of proceeds then as follows: 
a) Using (2.10) one gets 



< e 



Indeed, the r.h.s. is the i^-th Fourier coefficient of Q(ip) = rit)6efi'(')^ + i^Sv)', such a function is clearly 
analytic in the strip |Imi^| < p, and is there bounded by |tJ(',c)| < 1, so the estimate follows.^) 
b) From (2.13), by using the relation niy = h — I, one gets 

vee vee 



3) In the non analytic case, one would simply get one, in place of e P\^\ ^ at the r.h.s. of the inequality. In such special problem 
the improvement of the estimates due to analyticity is clearly irrelevant. 
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all integrals are then extended'*-' to Rett, G JR, and (2.14) is used. 

c) Finally, the sum over a introduces in the estimate a factor nt)(if + l)i 1^^^ than 2^^'^ , while the sum over 6 
gives a factor Mh = cardinality of 0^ ; one easily recognizes that Mh is precisely the quantity defined and 
estimated in Lemma 2, so one has in particular Mh < 4''"*. 

Putting everything together, one obtains 

\ah,:.\<Wi8CW)''-H-'">+^e-\''\''<^-^\ a = ^, q , (3.15) 

as claimed. 

Item iv) of the proposition is a particular case of (3.13), (3.14). 
3.3. The cancellation mechanism, and proof of point Hi). 

We show here that in the Hamiltonian case, due to a very special cancellation mechanism, the average 
energy rjhfi, which according to the above general estimate (3.15) does not contain any exponentially small 
factor, is instead proportional to e~^^'^'--^~^'> (it is interesting to note that a similar property does not hold 
for ^ft,o)- 

Let us consider again (3.13). The generic integral on ty ranges there from — cxd to tyi, v' being the vertex 
which immediately precedes v in the tree; the only exception is the integral associated with the root vertex 
(i.e. the vertex at the end of the root), which extends from — cxd to +cxd. We then use the trivial identity 

' 1 /■+°° 1 /■+°° 

.dty=- .dty+- . S{ty. - ty) dty , (3.16) 

) ^ J — CO ^ J — CO 

where S is the sign function, and distinguish the two terms of this decomposition by a label At, = for the 
former. At, = 1 for the latter. We then put a label At, , as a further decoration of the tree, on the branch 
leading to vertex v (the root has necessarily A = 0), and add in the formula for the value of the tree a primed 
sum over A (a sum over all v but the root); so we get 

V(0,a,n) = J2'v(0,a,n,\) , a^,. = ^ ^' ^ ^V((?, a, n. A) , (3.17) 



A 



with 



/.oo /.oo 
^ n -/ — CO J — CC ^ a 



Aion-i 



]^[\yS{ty -ty,) + {l-\y)] c'^^'-dt 



(3.18) 



It is very relevant to note that a label At, = interrupts the chain of integrations, and the value of the tree 
factorizes, with one factor for each cluster of vertices internally connected by branches with At, = 1. For 
an intuitive picture we think that branches with At, = are fragile and break, thus disconnecting the tree 
9 into I < h — 1 clusters Co, • • •,C;_i; let Co be the cluster connected to the root. As before, we can now 
raise the integration paths to Imtt; = ±r(l — b): the sign function S is not analytic, but this does not cause 
any problem if the sign of Imtt; is the same for all vertices of the same cluster; the sign of Imtt; can instead 
be different for different clusters, and we can profit of this new freedom. Indeed, for each cluster Ci let us 
denote 

Vi = ^ (3.19) 

f ec, 

{vq is called free momentum in [Gl], [G2]) so that 'Yl!i=o^i = = 0; by appropriately choosing the sign of 
Imty , we then get 

l-l 



A^ = ^|^^,| (3.20) 



Here the estimate is not optimal: by symmetrization of the integration domain, one could gain a combinatorial factor 
]^^(it'0~^' the final improvement on the constants A and B would be not so relevant. 
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Let us now introduce the decomposition 



"ft,o = a'ft.o + ) a=Lri, (3-21) 

where a'^ g collects the contribution of decorated trees with jyg 7^ 0, while a'^' g collects the contribution of 
trees with i>q = 0. If = but i>o 7^ 0, then clearly N > 2. Proceeding just as we did above to prove (3.15), 
one immediately obtains 

la'ft.ol < AB''-H-'"'+^e-^'^^^^-^^ (3.22) 

(the only new fact is that each sum over At, =0,1 gives now an extra factor 2, which however is compensated 
by the factor ^ in (3.18)). 

We shall now restrict ourselves to a = ry, and use the Hamiltonian character of the problem. Our aim is 
to show that one has exactly ry^' g = 0; in virtue of (3.22) this is clearly enough to prove point iii) of the 
proposition. To make evident the compensation, we must enumerate the trees differently. We now think 
of a tree with h vertices vq, ■ ■ ■ , Vh-i (before decorations are added) simply as a partially ordered set of 
h symbols, and count exactly one tree for each "tree-like" partial ordering one can introduce in the set 
{0, ...,/} — 1}. Tree-like will mean that the following particular rules are satisfied: 

i) one vertex (the root vertex) precedes all other vertices; 

ii) if both v' and v" precede v, then either v' precedes v" or v" precedes v' . 

Note that a branch drawn from v' to v indicates that v' immediately precedes v (any other vertex preceeding 
V also precedes v' ). 

The correspondence between old and new enumeration is seen after some meditation. For example, the 
three-vertices trees shown in figure 4a are now counted separately, while the two four-vertices trees exhibited 
in figure 4b are counted only once (they represent the same partially ordered set, simply drawn differently). 
The new set of trees with h vertices will be called 0^. Concerning decoration, the only difference is that the 
labels a = ^ or ry on the internal branches are now completely free (there is no more top-to-bottom order); 
the n and A decorations are left unchanged. 





(a) (b) 

Figure 4- Illustrating the new counting of trees as partially ordered sets: (a) different trees, to be 
counted separately; (b) different drawings of the same tree, to be counted only once. 

As a result of the new way of counting (for a general reference on the combinatorial problems concerning 
various trees, see [HP]), one gets for the value of trees expressions almost identical to the previous ones, in 
particular to (3.17) and (3.18), with precisely the following differences: (a) the Taylor coefficients f^J'-'" 
loose the combinatorial factor 1/ jy\(my — jt,)!, namely they are replaced by the derivatives 

la^ . ^3 23) 
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(b) an overall external factor 1/(1 + my)] = l/h\ is added. So one has 

1 
h 



^^••^ = t^J2T! J2T! ^(^' A) (3.24) 



with 



vee -^-^ -^-^ «6e ^ ' . (3.25) 

r^[XyS{ty -ty.) + {l- Xy)] e'^^'-dt 

We now focus the attention on the set of trees contributing to ry^' g, and divide it into equivalence classes; our 
purpose will be to show that one has exact compensation within each class. To define classes, we introduce 
the operation of "moving the root" of a tree (the idea is taken from a result of Chierchia taken from [Gl]: 
see the cancellation mechanism IV in §6 of [Gl] and the related lemma). First of all, let us note that in any 
tree the path from the root vertex, say vq, to any other vertex, say vi, is uniquely defined. Moving the root 
from vq to vi simply means modifying the order relation (i.e., replacing the tree inside by a different one) 
just by inverting the orientation of branches along the path. The idea, actually very intuitive, is illustrated 
in figure 5 (the change in the decoration will be explained in a moment). Note that the operation is well 
defined since, for a given set of (non oriented) branches, the choice of the root vertex (of the first element of 
the tree-like partially ordered set) uniquely determines the orientation (the partial ordering). 



t C 1 C 2 r] 1 C 2 




Figure 5. Illustrating the operation of moving the root in a tree. The root moves along the 
branch joining vertices and 1; correspondingly the branch changes orientation, and its label 
^ turns into rj. The second and the third graphs are simply different drawings of the same 
tree. 

We now say that two decorated trees {0,a,n,X) and (0' ,a' ,n' are equivalent, if: 

a) 6' is obtained from 6 by moving the root inside the root cluster Co- 

b) The decoration a' differs from a along the path: namely any label ^ is there converted into rj, and conversely 
(see figure 5); in compact notation, say that the label aj of branch h in the path is converted into its 
Hamiltonian conjugated, denoted aj. 

c) All other decorations are left unchanged, so n = n' and A = A'. 

To illustrate the compensation mechanism, it is convenient to introduce the following notation: 
Bo = the set of branches of Co 

v'^, v'-l = the ends (tail and tip respectively) of branch & G 5o- 

ko = the number of ry-type branches in Bo- 

Ioj Ho — restrictions oft and n to the cluster Co- 

d^a ^ = the derivative with respect to a = ^ or ry, which however operates only on /'-*'-'. 
We can write the value V oi a, decorated tree in the factorized form 

V{6 , a, n, A) = v''°\e, a, n) , a, n. A) (3.26) 

where V^'^^ is the contribution of the cluster Co, while V^'^^ collects the contributions of the remaining 
clusters; V^'^^ is clearly constant within the equivalence class. Using (finally!) the Hamiltonian character of 
the problem, namely definition (2.9) for and , one recognizes that for each branch h one has exactly one 

pair of derivatives d^a^^ and ±d^-'^\ with the minus sign if aj = rj; precisely, one has 

/. + 00 /. + 00 

!/(")((?, a, n) = f n 9nA / • • • / U{e,a,rv, i,) e'^^^ At^ (3.27) 
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with 



t^(^?,«,ffiio)=(-i)'° [ n n^«^^^4?^ Uf^^^i^^'^") (3-28) 



2 

heBo heBo veCo 



where is /, or a convenient derivative of / in case there are branches with A = coming out of v. Let 
us stress that dly "^ acts only on /'-*'-'. The outermost derivative 9^*'°-' refers to the root vertex, here denoted 
vq (accordingly to (2.9), the derivative is with respect to ^, since the label of the root is here rj). 
We can now see that 

U(0, a, n-t_o) = 5^"°^ U{e, a, n;to) , (3.29) 

with U independent of the tree within the equivalence class. Indeed, whenever one moves the root from one 
vertex to a nearby one (it is clearly enough to consider this case), a branch h of Co is inverted: correspondingly 
iij e t;" (tail and tip of h) interchange, but on the same time aj is replaced by aj (and conversely), so the 

product da^ d—* is left invariant. The sign function S relative to the inverted branch changes its sign, but 
the number ko of ry-type branches changes exactly by one unity, so U does not change. 
The conclusion now follows: by taking the sum over all trees in the same equivalence class, one clearly gets 

( E ^t ^) ^(^' a> ffi io) = + • • • + ^) l^iS, a ffi h) (3.30) 

veCo ^ ^° 

where each derivative ^j- is intended to act only on inside U , and not on sign functions. One can 
finally integrate by parts: for the assumptions on / at infinity the finite term vanishes; the derivatives of 
the sign functions give no contribution (more precisely, one easily sees that each of them gives two opposite 
contributions); and when the derivatives are transferred to the exponential, one gets in front an overall factor 
'^vec = = 0. This proves the exact compensation inside a class, so that ry^' g = 0, and concludes the 
proof of the Proposition. 



4. Concluding remarks. 

4.1. I-dependent, non factonzed perturbations. 

The Hamiltonian (2.1) that we have treated, has a rather special form; in particular, the coupling term 
eg{ip)Vi{x) is factorized and independent of I. Such Hamiltonian, properly speaking, does not represent 
the motion of an oscillator in an external potential, rather a point mass which is periodically perturbed 
with fixed frequency to. The consequence of our choice was that, once the coordinates ^, rj were introduced, 
the equations of motion assumed the special form (2.8), namely they reduced to only one pair of (time 
dependent) equations for ^ and rj. For this reason we could write a perturbative scheme to compute ^(t) and 
ri{t) at all orders in e, treating separately (and trivially) the motion of Lp, and completely ignoring I. 

However, the perturbative scheme we developed adapts very easily to Hamiltonians of the form 

A'(e, ry, ^, 7) = + ry + ef{i, ry, ^, I) , (4.1) 
or more generally, to equations of motion of the form 

with, in the Hamiltonian case, = |^ , = , = |^ , fj = 

We indicate briefiy how to proceed in this more general situation, and to this end, we introduce the following 
compact notation: 

x = {tv,v,i) 

a = a label running on the set A = rj, (f, 1} 
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i = {jeJv'jv'ji) e ff^; lil = Eai«; and j! = j^! • • -jj! 
t={k^ eIN+:aeA, 0<i< i„}, and |i| = ^ fc" 

X(t) = Xo(t) + ^£''X,(t) 
ft 

XL = . . . pi 
^" ~ i! 5X1 • 

The recursive scheme that one gets, for a perturbative representation of the motions satisfying the asymptotic 
conditions 

C(t)-t^O, v(i)^v\ ^(t)-iot^^\ I(t)^I" (4.3) 

for t —CO, is 



«ftW = / mnn' 

j:l<|j|<ft-l fc:|fc|=ft-l l<ii<ji l<ii<ji 



i 



Note that \j\ replaces here m. 

Such a scheme gives rise to tree expansions quite similar to the previous ones; the only difference is that now 
the labels a on root and branches run over the whole set A, and not only on ^ and rj. One then introduces, 
as before, the Fourier expansion of the perturbation, here 

/^(e,ry,^,7)= E^«.^(€,'?,^)e'''^ (4.5) 

and collects terms which contribute to the same Fourier component a^^^, a ^ A; this means introducing 
sums in the perturbative expansion, and adding decorations to the tree, exactly as before. The estimates 
then run quite smoothly, and one finds exponential estimates, i.e.: 

\ah,^\ < AB''-H-'"'+^ p-\^\u>r(i-s) ^4 g^, 

precisely as before. An assumption generalizing (2.10), (2.13) is 

|j^^,^(s + i(T,??°,r)| < Clil(5-lil-™° w;(s)e-^l''l for |(t| < r(l - ,5) (4.7) 

In the Hamiltonian case, and only for a = rj or a = I (as is obvious, Arj and A7 are not independent, since 
for energy conservation one has asymptotically Arj = — wA7), one has then the cancellation mechanism, 
leading to the exponential laws 

\ah,o\ < AB''-H-'"'+^ e-'^'^^^-"'^ , a = 7], I (4.8) 

The method is exactly the same: one uses the identity (2.23), and introduces clusters of vertices as before; 
everything repeats then precisely as in the previous case, with the only obvious generalization that, when 
the root is moved along a path in the cluster Co, labels ip and I also interchange. The conclusions, essentially 
expression (3.17) for rj and a very similar one for I, follow. In a word: the case of Hamiltonian (4.1) is 
practically identical to the simpler case we treated in detail; only the notation gets complicated. 



4.2. Reversible systems. 

The cancellation mechanism leading to the special exponential law for rjhfi was based in an essential way 
on the Hamiltonian character of the equations of motion. It is worthwhile to remark that an alternative 
assumption leading to the same result is a symmetry in the system, namely that f^iCrj) and frjiCf]) are 
respectively even and odd in ^, and gi^p) is even in ip. The time reversal symmetry of (1.1), t.e. the 
transformation i : (px,Py,x,y) (—px,—Py,x,y), becomes in the ^,ri,ip,I variables i' : (^,ri,ip,I) 
(— ^, rj, —ip, I): hence the above parity properties imply time reversal invariance of the equations of motion 
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(i.e. they imply that i' St = S-ti' , ii t ^ St is the flow generated by the equations of motion). Hence 
reversibility can replace the Hamiltonian character. 

In such a non-Hamiltonian but reversible cases the proof that Tf\fi = proceeds as follows. First of all, it 
is not necessary to change the enumeration of trees (i.e., redeflning trees as ordered sets and so on). Each 
tree can be described as a core, namely the root cluster Co fully composed of branches with A = 1, carrying 
several (possibly none) "sub-trees" connected to the core by branches with A = (the roots of the sub-trees). 
Instead of moving the root, here we move sub-trees, more precisely sub-trees with ^-type root, if any. Let 
us sketch the procedure. Assume for a moment that at least one ^-type sub-tree exists^ and has multiplicity 
one, namely it is different from all other sub-trees, and consider the "pruned tree" 9 one obtains by just 
removing this sub-tree. The value V{9,a,n,X) is obtained by integrating the functions (see (3.18)): 



) \[\S{U + e'^^'- . (4.9) 



ve8 

By adding the sub-tree to any vertex of Co, in all jo + 1 possible ways, one gets 



[Yifar'^iiv,^")] (Jv + l)fZ'^'-''+HU,vl lll[KS(U-U.) + (l-K)]e''^^^ (4.10) 



But clearly, one has 



ih + i)fa;^''''+' = ^fa;'' , (4.11) 



so that, also taking the sum for G Co, one gets 

+ • • • + ^) n fa:'''iiv,vi kKS{t. - 1^,) + (1 - A,: 



(4.12) 



where the derivative acts only on f^J'-'" (and not on sign functions containing t^). Integrating by parts 
as before, one easily gets a'^' ^ = 0, a = ^, ry. 

In the case there are ^-type subtrees, but all of them have multiplicity greater than one, one can still move 
any of them, say of multiplicity k; simple combinatorial arguments then show that one gets 

/ d d \ 

in place of ^ + . . . + ^jjj— in the above expression, and the conclusion does not change. 

The operation of "moving a sub-tree", here introduced, closely resembles the operation of "moving a fruit" 
in [Gl] (see mechanism II in §6 of [Gl]). 

So far the symmetry of / and g has not been exploited;®) we use it now, to take care of trees without any 
^-type subtree. One easily checks that, if the root is of rj-iype, then the quantity to be integrated is exactly 
antisymmetric for the exchange to ~^ ~io, and simultaneously no ~^ "Ho, in root cluster. Indeed: 

i) The number of S functions which change sign is /lo — 1. 

ii) The number of derivatives with respect to ^, which obviously change the parity of /„ functions, is XlugCo • 

iii) The number of odd functions (including the one associated with the root) is 1 + '^y^Cgi^^ ~ jv)- 

The overall parity is thus hg + Xlugc "^f ~ ^/lo — 1, namely odd. On the other hand, the scalar product 
n ■ t is invariant, while, by the assumed parity of g, one has g„ = g-n- Since all integrals run from — cxd 
to cxD, one concludes that trees with opposite decorations n and —n have opposite value, thus they exactly 
compensate when one computes ry^' g. The latter cancellation mechanism is close to the "parity cancellation" 
(see mechanism I in §6 of [Gl]). 

In essence: the above analysis puts in evidence the terms which are possibly large, among those contributing 
to Arj, namely those associated to trees with = 0, i^o = 0, and no ^-like branches issuing from Co. Such 
terms cannot be proven to be small on the basis of purely perturbative arguments (they would not be small. 



^) Hence the above argument has general vahdity, and it could also be used to (slightly) improve constants in the estimates of 
Section 3. 
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for example, for a problem with friction): they compensate exactly only in presence of a symmetry, namely 
either the Hamiltonian character or the reversibility of the equations of motion. 

4.3. The problem of the energy exchange as a problem of split separatrices. 

Here we provide a (rather informal) geometric picture, which can be useful to interpret our results; in 
particular, we establish the link between the problem of the energy exchange, which was discussed in the 
previous sections, and the idea of separatrices splitting. 

To this purpose, let us consider again the Hamiltonian K given in (2.6), and look at the three-dimensional 
energy surface Te corresponding to some fixed value of the total energy E. Good coordinates on F^; are 
(^,ri,ip), but for an easier insight it is better to look at ip° = ip — uit in place of ip. We fix the attention 
to the set of trajectories in Te, such that ^(t) — t ^ asymptotically for t — cxd.®) One has then one 
trajectory (one curve in Te) for each choice of rj" = lim(^_oo f](t) and of ip° . By varying ip° at fixed rj" , one 
gets a two-dimensional regular submanifold S~„ ^ of F^; (in fact, an analytic surface), whose equation can 
be written in the form rj = ry~„ ^(^,1^°). Clearly S~„ ^ coincides, for e = 0, with the plane (the cylinder, if 
one recalls that ip is an angle) of constant rj = rj" . For e ^ 0, S~o ^ is (by definition) fiat asymptotically 
for ^ —CO, while, according to Proposition 1 (Section 2), it assumes a sinusoidal shape, of exponentially 
small amplitude, asymptotically for ^ +oo; in the middle, namely in correspondence with the collision, 
the oscillations are instead much larger, namely of order e. The behavior of S~o ^ is represented in figure 6 
(gray surface). 




Figure 6. Illustrating the behavior of ^ (white) and S^o ^ (grey). 



Exactly in the same way, one can look at trajectories such that ry ^ ry° in the opposite limit t +cxd. The 
corresponding curves form a second surface ^, described by an equation of the form rj = ry^„ ^(^,1^°), 
which is instead fiat for ^ +00, and assumes a sinusoidal shape (of exponentially small amplitude) for 
^ —CO, with oscillations of amplitude 0(e) in the middle (figure 6, white surface). We claim that the two 
surfaces intersect each other as is symbolically shown in figure 6, more precisely: 

i. By taking a section ip° = const., one gets two curves a^o and a~o which either coincide (this 
happens for those special values of ip° , such that the asymptotic energy exchange Arj vanishes; there are two 
of them for any rj"), or never intersect. 

Indeed, should they intersect once, they would intersect infinitely many times, and this is in confiict with 
the existence of well defined limits of both curves in both directions. 

ii. Although rj^o ^ and ry~o ^ are in a sense large (namely of order e) during the interaction, their difii"erence 
at fixed ^ and ip° is always exponentially small: 

^p,ALvn-^;o^Ai,^n = Oie-n (4.14) 
so, in particular, a generic section (f^ = const, looks like figure 7. 



^) The choice is made for coherence with the previous notations; actuaUy, one could equaUy weU make the more symmetric choice 

C(o) = o. 
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Figure 7. Possible profiles of a^o ^ and a^o ^ • 

This follows from volume preservation of Hamiltonian dynamics: indeed, in place of K , one can equivalently 
consider, for the evolution of ^ and rj, the non-autonomous Hamiltonian 

K(i, rj, t; ^") = v + eg{^" + ujt)f{^, rj) (4.15) 

depending parametrically on ip° . The area d^dry is clearly preserved by the Hamiltonian dynamics; on the 
other hand, d^ is almost preserved, as follows from = ^ dt = (I + 0(e)) dt. So dry is almost preserved, 
too, and any vertical segment in figure 6 — in particular, the vertical separation of a^o^^(ip°) at any ^ — is 
in turn preserved by the dynamics, up to a factor 1 + 0(e). 

iii. In any section ^ = const, one obtains split curves a^o ^, which intersect transversally, but for any ^ their 
separation is exponentially small (see figure 8). 



(a) (b) (c) 

Figure S. Possible profiles of (7~o ^ (solid line) and (T^o ^ (dashed line) at ^ = — oo 
(a), during the collision (b), and at ^ = +00 (c). 

This is an obvious consequence of the previous analysis, which however shows an important point: the 
problem of the energy exchange in our system can be thought of as a problem of separatrices splitting. This 
is relevant, in our opinion, for the connection with [G1,G3] (as a matter of fact, this was our first approach 
to the problem: only later we recognized that the direct computation of the energy exchange is simpler). 

The special form of the perturbation in the Hamiltonian (2.6) is not relevant, and essentially identical 
conclusions can be drawn for the more general case considered in (4.1). 

4.4. Systems of several oscillators. 

An interesting question is whether the above analysis and results can be extended to systems with more 
than one fast degrees of freedom. Work is in progress in this direction. We anticipate here just one result, 
concerning a natural extension of Hamiltonian (2.6), namely 

K(C,n,^,L)=^-L+n + ef(C,n)9(^) , (4.16) 

with T], Q IR, ip = (ipi, . . . , ipi) G it', 7 = (7i, . . . , /;) G iR'. Such a system represents a point mass which 
is perturbed in a quasi-periodic way, with / frequencies. 

We keep the assumptions (2.13), (2.14) for /, and trivially generalize (2.10) to |(j'(',c)| < 1 for |Imi^j| < p, 
j = I, . . . ,1. Proposition 1 can be extended as follows: 

Proposition 2. Consider the Hamiltonian (4.16): 1st f and g satisfy the above assumptions, and consider 
motions satisfying the asymptotic conditions (2.15). Then: 

i) For small e the Taylor series (2.17) converge, precisely for any real t one has 

\an{t)\ < AB''-^ « = (4.17) 
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ii) For any v_ G 'E , the Fourier component ah^v of ah{-\-(X}), a = rj, satisfies the exponential estimate 

\ah,,\ < A5''-i,5-''«+ie-^l^'^l(i-'^)-^l^l . (4.18) 

Hi) For !E' such that v_ - UJ_ = 0, one has 

\r^n,u\ < A5''-i,5-''«+ie-2^l^'^^l(i-'^)-2Hii'l . (4.19) 

with some i/ such that il' • w 7^ 0. 

iv) At first order in e one has the explicit expression 

/oo 
f{t,r,»)e'\"-^\' At . (4.20) 
- 00 



Since v_-uj_ can be arbitrarily small for large v, the therm — at exponent in (4.18) is now very important. 
The exponential estimate (4.19) is the analog of (2.22), and is based on cancellations. Expression (4.20) 
represents the JLT approximation in this problem. Estimates of the form (4.18) are also found in [Gl] (see 
(8.1)) for separatrix splitting problems (they were, however, derived with a different technique relaying on 
the classical proofs of the KAM theorem), and in [Ge], where they are derived by using a tree representation 
of the ah^v See also [Si], where a first order analysis is performed in connection with the related problem of 
separatrices splitting. 

For a detailed analysis of the JLT approximation in systems with several fast variables, including a careful 
comparison with numerical results (quite nice agreement is found), see [BCF]. 

4.5. Improving the bounds. 

The presence of the factor {1 — 6) in the exponential bounds (2.21) is clearly a sign that some of the estimates 
can be improved. Of course one needs more information about the nature of the singularity at Imt = ±i. If 
we keep the assumption that Vq = and Vi = then we see that have a pole of order mo = 3 in 

the t variable at ±ir. 

This means that the integrals on the variables can be shifted to larger imaginary parts than r, leaving, 
besides the integrals at Imt^ = r (1 + (5), a few integrals over segments ±[ir, ir (1 + b)\ (due to logarithmic 
singularities generated from the poles by the integrations) plus some residues from the singularities at ±ir. 
Since there are h integrals the residues, which provide the most singular contributions, will contribute 
quantities bounded by: 

(C^|n,|w)(™°-i)''e-l''l'^^ (4.21) 

V 

Hence if g{ f ) is a trigonometric polynomial of degree N our estimates can be improved, with some extra 
work, to: 

/i''(w#C)2''e-'^^(l''l+2'^°-') (4.22) 

This shows, by combining the latter bound with the previous and by optimizing the choices, that in the 
above cases the condition for having domination of the lowest order is essentially (up to logarithms) eoj^ <^ 1 
(instead of the su)^ <^ 1) discussed in §1. 
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The Appendix 



As already remarked in the Introduction, the use of the ^, rj variables simplifies the perturbative scheme, but 
is not strictly necessary. We give here a few indications on how to proceed with the original p, q variables; this 
may be useful on the one hand for a comparison with [Gl] and [G2], on the other hand because, although the 
scheme gets more complicated for the presence of nontrivial Wronskians and kernels in integral expressions, 
some points, like the analyticity requirements on the perturbing function, are perhaps more immediately 
understood. 

The system we aim to study is Hamiltonian (2.1), with exponential potentials. However, to make more 
transparent the perturbative construction, and emphasize in particular the analyticity assumptions we need, 
we write the equations of motion in the more general form 

i = p, p=<^{x) + sg{(p)f{x), ip = uj, i = ■^{'f)V {x) ; 

in the problem at hand one has clearly = f{x) = V{x) = e~^, and '){x) = —g'(ip). We look for motions 

such that, asymptotically for t -co, 



p{t) ^ p{-oo) = -^/2r], x{t) - p{-oo)t ^ 0, Lp{i) - Loi Lp\ I{t) ^ I{-oo) = 1° ; 

we aim to compute A7 = 7(+cxd) — 7(— cxd). 
For £ = 0, and = e~^, one has the unperturbed motion (see (2.5)) 

xoit) = 2\og[r]-\cosh^/^tf], poit) = ^/2^tanh^/^t, Lpo{i) = Lp° + Loi , Io{t) = T , 

and the associated Wronskian, restricted to the first two components, is 

y^U)=f^^^i^) ^^pi^)\ ^ ( 1-yfttanhyft y|tanhyfA 

\wp,it) Wppii)) \_tanhyft- §[coshyft]-2t [cosYi^t]-'' ) ' 

For £ we write 

a;(0 = a;o(0 + ^£''«ft(0> =Po(0 + ^£Vft(0> /(t) = r + ^£''7ft(t) , 

ft>0 ft>0 ft>0 

and all series will converge provided, in the series for x{t), one extracts at each order h a common factor j3{t) 
growing linearly with \t\: say, to fix the ideas, (3{t) = \/l -\-t'^ (such a linear growth is due to the element 
Wxx in the wronskian above). 

We then expand $, / and V in Taylor series around the unperturbed motion xi:j{t): 

^x{t)) = Hx,{t)) + Y.'^'"^Mm<i)-Mt)r , with <i>™ = i^|^, 

m> 1 

and similarly for / and V . One must require that such series converge, although, as remarked above, 
x{t) — xi:j{t) diverges linearly in t. The simplest situation is that of integer functions, like our exponentials: but 
functions with singularities are also welcome, provided the radius of convergence of their Taylor expansions 
depends on the point of expansion xq, and tends to infinity for xq oo. A very natural natural assumption, 
although not the most general one, for the scattering problem we are dealing with, is that smgulanttes of 
all functions have hounded real part. As in the previous approach, we also need, of course, that all functions 
decay to zero (faster than x~^) for x ^ oo. 

By substituting the series expansions in the equations of motion, we get, for any h> Q, 

Xh(t)= I KAt,t')Fh(t')dt', Ph(t)= I Kp(t,t')Fh(t')dt', Ih(t)= I F(t')dt' , 

J — CO J — CO J — CO 

where K^, Kp are antisymmetric kernels, 

Kx(t,t') = -Wxx(t)wxp(t') + Wxp(t)wxx(t'), Kp(t,t') = -Wpx(t)wpp(t') + Wpp(t)wpx(t'), 
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and Fh, Fh are given by 



2<m<h fel,---,fe„^>l 
\k\=h 



l<m<h-l fel,---,fe„^>l 
\k\=h-i 

Fhit) = -ji^° +iut) J2 ^"(^oW) J2 ^kAi)---^kJt) . 

l<m<h-l f.i,^^^,f.™>i 

\k\=h-l 

For h = I such expressions should be red 

Fi(t) = g(^" + iot)f(xo(t)) , Fi(t) = + iot)V(xo(t)) . 

^From these relations, one naturally works out tree expansions rather similar to the previous ones; the 
only difference is that now one has only one type of branches, namely branches "of type x" , but two types 
of vertices, depending on the choice of the term with $"* or with /"* in the expansion for Fh- We shall 
distinguish the two possibilities by associating to any vertex v a label jy = or jy = 1, respectively. At 
variance with the previous case, only vertices with jy = 1 now contribute to the order in e, more precisely 
one has h = jy. An example of a tree for x, more precisely a tree contributing to xg, is reported in figure 
9; open circles and dots correspond there respectively to jy = and jy = 1. Formal rules are that (i) any 
end vertex v has jy = 1; (ii) niy > 2 if jt, = 0; (iii) trees for Ih have j = 1 at the root vertex. Trees for Xh, 
Ph and Ih differ only for the root. 




Figure 9. A tree with mo = 2, mi = 2, m,2 = 3, ms = 2, m4 = 2, . . .; open circles and dots 
correspond respectively to =0 and = 1 . The tree represents a contribution to ^rg (t) . 



As in Section 3.2, each tree with m vertices represents a m-tuple integration over variables dto • • ■'^tm-i, 
with integrals nested accordingly to the order of the tree; the function to be integrated is, here too, a product 
of one factor for each vertex: precisely, any vertex v but the root vertex contributes with a kernel K^iiv' ,iv), 
v' denoting the vertex immediately before v, while the root vertex vq has kernel Kx{i,ivo), or Kp(t,tyg), 
or no kernel at all, respectively for Xh, Ph or Ih- Moreover, any vertex with jy = contributes with a 
factor (^"^^ (xo(ty)), while, for jy = 1, the factor is g(ip° + uity)/"^^ (xo(ty)) or, for the root vertex of /^(t), 

+iOty,)V'^^o(xo(ty,)). 

Everything proceeds then smoothly in the same way as with the energy-time variables, namely: 

a) One makes the position: = sup^^jj^ \xh(t)\/ l3(t); proceeding as in Section 3.1, it is not difficult to find, 
here too, a recursive relation of the form 

h-l 

Ch < const. ^ Ci Ch-i , 
1=1 
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which assures convergence. 

b) Exactly as in Section 3.2, one expands g and 7 in Fourier series, and further decorates trees by indicating, at 
each vertex v with jj, = 1, a Fourier index ny G !E. Vertices with jy = have, by definition, ny = 0, so the 
Fourier index, or momentum, of the whole tree is = ny. As in Section 3.2, one raises simultaneously 
all integration paths, getting, for Xh, Ph and Ih, and for nonvanishing momentum, exponential laws of the 
form (3.15). 

c) Finally, to prove that the average (momentum = 0) of Ih follows the special exponential law (3.22), one 
produces relevant cancellations. The procedure closely follows Section 3.3, so we shall not describe it in 
detail. Basically, here too one introduces the decomposition (3.16), which leads to clusters, and proves 
exact cancellations among trees with vanishing free momentum (i.e., with Fourier component no of the 
root cluster): the cancellation mechanism is, in this case, exactly the same as the one pointed out by L. 
Chierchia and exposed in [Gl]. The equivalence classes giving exact compensations are obtained, as in the 
previous approach, by moving the root from one vertex with j = 1 to a similar one inside the root cluster; 
by the way, just here one uses the Hamiltonian character of the problem, namely the fact that 7 = —g' 
and / = —V'. A crucial point is that the product K^iiv' ,iv)S{tyt — ty) is symmetric (in other words: the 
antisymmetric factor plays the same role as the change ry ^ ^ on branches in the path, in the previous 
approach). 



References 

[ACKR] Amick C, Ching E.S.C., Kadanoff L.P., and Rom-Kedar V.: Beyond All Orders: Singular Pertur- 
bations m a Mapping J. Nonlinear Sci. 2, 9-67, 1992. 

[BCF] Benettin, G., Carati, A., and Fasso, F.: Adiahatic invariants for many -frequency systems: a careful 
comparison between perturbation theory and numerical results, paper in preparation. 

[BCS] Benettin, G., Carati, A., and Sempio, P: On the Landau-Teller approximation for the energy ex- 
changes with fast degrees of freedom, Journ. Stat. Phys. 73, 175-192 (1993). 

[BGGl] Benettin, G., Galgani, L., and Giorgilli, A.: The dynamical foundation of classical statistical 
mechanics and the Boltzmann-Jeans conjecture, in Seminar on Dynamical systems, edited by S. Kuksin, 
V.F. Lazutkin and J. Poshel (Birkhauser, Basel 1993). 

[BGG2] Benettin, G., Galgani, L., and Giorgilli, A.: Realization of Holonomic Constraints and Freezing of 
High-Frequency Degrees of Freedom m the Light of Classical Perturbation Theory. Part I, Comm. Math. 
Phys. 113, 87-103 (1987); Part II, Comm. Math. Phys. 121, 557-601 (1989). 

[Bel] Benettin, G.: On the Boltzmann-Jeans Conjecture for a Classical ''Freezing" of Fast Degrees of 
Freedom, proceedings of the 1991 Taormina Conference on Thermodynamics, edited by P.V. Giaquinta, 
M. Vicentini Missoni and F. Wanderling (Accademia Peloritana dei Pericolanti, Messina 1992). 

[Be2] Benettin, G.: Adiabatic invariants and time scales for energy sharing m models of classical gases. 
Proceedings of NATO conference on Hamiltonian systems, Toruri 1993, edited by J. Seimenis (Plenum Press, 
to appear). 

[BF] Benettin, G., and Fasso, F.: Classical ''freezing" of plane rotations: a proof of the Boltzmann-Jeans 
Conjecture, Journ. Stat. Phys. 63, 737 (1991). 

[Bol] Boltzmann, L.: On certain questions of the theory of gases. Nature 51, 413 (1895). 

[Bo2] Boltzmann, L.: Lectures on gas theory, translated by S.G. Brush, University of Cal. Press (1966); see 
especially Part II, Section 45, Comparison with experiments. 

[C] Carati, A.: // metodo di Jeans-Landau-Teller per lo studio degli mvarianti adiabatici, PhD Thesis, 
Padova 1994. 

[CBG] Carati, A., Benettin, G. and Galgani, L.: Towards a rigorous treatment of the Jeans Landau Teller 
method for the energy exchanges of harmonic oscillators, Commun. Math. Phys. 150, 321 (1992). 

[CFl] Chierchia, L., and Falcolini, C: A Direct Proof of a Theorem by Kolmogorov m Hamiltonian Systems, 
Annali della Scuola Normale Superiore di Pisa, 21, 541-593 (1994). 



29agosfol995 



22 



[CF2] Chierchia, L., and Falcolini, C: Compensations m small divisor problems, preprint, 77" Universita 
di Roma, august 1994, archived in mp-arc@ math, uiexas. edu, #94-270. 

[CG] Chierchia, L., and Gallavotti, G.: Drift and diffusion m phase space. Annales de I'Institut H. Poincare 
60, 1-144 (1994). 

[E] Eliasson L. H.: Hamilionian systems with linear normal form near an invariant torus, ed. G. Turchetti, 
Bologna Conference, 30/5 to 3/6 1988, World Scientific, 1989. And Generalization of an estimate of small 
divisors by Siegel, ed. E. Zehnder, P. Rabinowitz, book in honor of J. Moser, Academic press, 1990. But 
mainly: Absolutely convergent series expansions for quasi-periodic motions, report 2-88, Dept. of Math., 
University of Stockholm (1988). 

[Gl] Gallavotti, G.: Twistless KAM ton, quasi flat homoclinic intersections, and other cancellations m the 
perturbation series of certain completely mtegrable hamiltonian systems. A review. Reviews Math. Phys. 6, 
343-411 (1994). 

[G2] Gallavotti, G.: Perturbation Theory, in "Mathematical Physics towards the XXI century", p. 275-294, 
ed. R. Sen, A. Gersten, Ben Gurion University Press, Ber Sheva, 1994. 

[G3] Gallavotti, G.: Twistless KAM ton. Communications in Mathematical Physics, 164, 145-156, 1994. 

[Ge] Gentile, G.: A proof of existence of whiskered ton with quasi flat homoclinic intersections m a class of 
almost mtegrable systems, preprint Roma "La Sapienza" n. 996, february 1994, archived in mp-arc@ math, 
utexas. edu, ^^94-33, in print on Forum Mathematicum. See also Gentile, G.: Whiskered ton with prefixed 
frequencies and Lyapunov spectrum, preprint Roma "La Sapienza" n. 1015, april 1994. archived in mp-arc@ 
math, utexas. edu, #94-92. In print in Dynamics and Stability of Systems, 1995. The latter papers are 
complete and correct even though the results, obtained with a very close formalism and notations, have later 
been stated in [CF2]. 

[HP] Harary, F., and Palmer, E.: Graphical enumeration. Academic Press, 1973, New York. 

[JKP] Joye, A., Kunz, H., and Pfister, Ch.-Ed., Exponential decay and geometric aspects of transition prob- 
abilities m the adiabatic limit. Annals of Physics 208, 299-332 (1991). 

[Jl] Jeans, J.H.: On the vibrations set up m molecules by collisions, Phil. Mag. 6, 279 (1903). 

[J2] Jeans, J.H.: On the partition of energy between matter and Aether, Phil. Mag. 10, 91 (1905). 

[JPl] Joye, A., and Pfister, Ch.-Ed., Superadiabatic evolution and adiabatic transition probability between 
two nondegenerate levels isolated m the spectrum, J. Math. Phys. 34, 454-479 (1993). 

[JP2] Joye, A., and Pfister, C: Semiclassical aymptotics beyond all orders for simple scattering systems, 
SIAM Journal on Mathematical Analysis, 26, 944-977 (1995). 

[LT] Landau, L and Teller, E.: On the theory of sound dispersion, Physik. Z. Sowjetunion 10, 34 (1936); 
also in Collected Papers of L. D. Landau, edited by D.Ter Haar (Pergamon Press, Oxford 1965), page 147. 

[MN] Martin, Ph., and Nenciu, G.: Semi classical inelastic S -matrix for one- dimensional N -states systems, 
preprint ESI #89, april 6, 1994, Vienna, Pasteurgasse 6/7, A-1090 Wien, Austria, and appended references. 

[Nl] Neishtadt, A.I.: On the accuracy of conservation of the adiabatic invariant, Prikl. Mat. Mekh. 45:1, 
80-87 (1981) [J. Appl. Math. Mech. 45:1, 58-63 (1982)] 

[N2] Neishtadt, A.I.: The separation of motions m systems with rapidly rotating phase, Prikl. Matem. 
Mekan. 48:2, 197-204 (1984) [J. Appl. Math. Mech. 48:2, 133-139 (1984)] 

[PV] Percival, I., and Vivaldi, F.: Critical dynamics and diagrams, Physica, D33, 304- 313, 1988. 

[R] Rapp, D.: Complete classical theory of vibrational energy exchange, Journ. Chem. Phys. 32, 735 
(I960). 

[S] Sempio, P: Scambi di energia tra gradi di libertd carattenzzati da differenti scale temporali: rilassamento 
all'equilibrio m gas poliatomici e problemi di confinamento magnetico , PhD Thesis, Padova 1993. 

[Si] Simo, C: Averaging under fast quasipenodic forcing, in Hamiltonian mechanics, mtegrability and chaotic 
behavior, J. Seimenis editor, Nato ASI Series B 331 (Plenum Press, New York 1994). 



29agosfol995 



23 



